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ABSTRACT 

The existence of outflow in aceretion flows is eonfirmed by observations and magneto- 
hydrodynamics (MHD) simulations. In this paper, we study outflows of accretion flows 
in the presence of resistivity and toroidal magnetic field. The mechanism of energy 
dissipation in the flow is assumed to be the viscosity and the magnetic diffusivity due 
to turbulence in the accretion flow. It is also assumed that the magnetic diffusivity 
and the kinematic viscosity are not constant and vary by position and a-prescription 
is used for them. The influence of outflow emanating from accretion disc is considered 
as a sink for mass, angular momentum and energy. The self-similar method is used to 
solve the integrated equations that govern the behavior of the accretion flow in the 
presence of outflow. The solutions represent the disc which rotates faster and becomes 
cooler for stronger outflows. Moreover, by adding the magnetic diffusivity, the surface 
density and rotational velocity decrease, while the radial velocity and temperature 
increase. The study of present model with the magnitude of magnetic fleld implies 
that the disc rotates and accretes faster and becomes hotter, while the surface den- 
sity decreases. The disc thickness increases by adding the magnetic field or resistivity, 
while it becomes thinner for more losses of mass and energy due to the outflows. 
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1 INTRODUCTION 

Accretion disc is an important physical object in astro- 
physics. The standard model of thin accretion disc is widely 
regarded as a successful model for explaining the observa- 
tional features in active galactic nuclei (AGN) and X-ray 
binaries (Shakura & Sunyaev 1973). However, the standard 
disk model cannot explain the spectral energy distributions 
(SED) of many sources, such as Sgr A* . To understand such 
systems, the model of advection dominated accretion flow 
(ADAF) is introduced (Ichimaru 1977; Narayan & Yi 1994). 
In this accretion flows, the energy released due to dissipa- 
tion processes is retained in the fluid rather than being ra- 
diated away. The models of ADAF place an intermediate 
position between the spherically symmetric accretion flow 
of non-rotating fluid (Bondi 1952) and the cool, thin disc of 
classical accretion disc theory (e. g. Pringle 1981). 

The observational evidence of accretion flows imply that 
outflow is an important property in these systems. For exam- 
ple, comparison of the accretion luminosities of neutron stars 
and white dwarfs in quiescence with similar binary com- 
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panions implies that the accretion rate onto white dwarfs 
is larger by three orders of magnitude than that on the 
surface of neutron star (Loeb et al. 2001). This indicates 
significant outflows in accretion flows. Moreover, outflows 
seem to be common in the nuclei of galaxies. Marrone et 
al. (2006) suggested the accretion rate of Sgr A* at small 
radii, much smaller than the Bondi radius, must be low, be- 
low 10"^MQj/r"\ While, Baganoff et al. (2003) estimated 
the hot plasma surrounding Sgr A* should supply the accre- 
tion rate of 10~^ Mqut"^ at the Bondi radius. This signifi- 
cant difference between the inner and outer mass accretion 
rates indicates mass loss from the accretion flow due to the 
outflow (Kawabata & Mineshige 2009). Thus, the accretion 
discs in the presence of outflow have been studied by several 
authors (Knigge 1999; Fukue 2002, 2004; Shadmehri 2008; 
Xie & Yuan 2008; Kawabata & Mineshige 2009; Bu et al. 
2009; Li & Cao 2009; Abbassi et al. 2010). 

Knigge (1999) derived the radial distribution of dissi- 
pation rate and effective temperature across a Keplerian, 
steady-state, mass losing accretion disc, using a simple, 
parametric approach that is sufficiently general to be appli- 
cable to many types of dynamical disc-wind models. Fukue 
(2002) examined a hydrodynamical wind, which emanates 
from an accretion disc and is driven by thermal and ra- 
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diation pressure, under an one-dimensional approximation 
along supposed streamlines. In another study, Fukue (2004) 
studied a supercritical accretion regime, where the mass ac- 
cretion rate was regulated just at the critical rate with the 
help of wind mass-loss. He derived a critical radius that out- 
side of it, the disc is in a radiation-pressure dominated stan- 
dard state, while inside this radius the disc is in a critical 
state, where the excess mass is expelled by wind, and the 
accretion rate is kept being just at the critical rate at any 
radius inside of critical radius. Shadmehri (2008) studied the 
effects of thermal conduction and outflows on ADAFs. He 
found that in comparison to accretion flows without winds, 
the disc which rotates faster and becomes cooler because of 
the angular momentum and energy flux which are taking 
away by the winds. Xie & Yuan (2008), based on one-and- 
a-half-dimensional description of the accretion flow, consid- 
ered the interchange of mass, radial and azimuthal momen- 
tum, and the energy between the outflow and inflow. Bu 
et al. (2009) presented the self-similar solutions for ADAFs 
with outflows and ordered magnetic flelds. They assumed 
the magnetic field has a strong toroidal component and a 
vertical component in addition to a stochastic component. 
They found that the dynamical properties of ADAFs can 
be significantly changed in the presence of ordered magnetic 
fields and outflows. Abbassi et al. (2010) examined the ef- 
fects of a hydrodynamical wind on ADAFs in the presence 
of a toroidal magnetic field under a self-similar treatment. 
Their results implied that in the presence of the wind, the 
disc temperature decreases due to energy fiux, which is taken 
away by winds and the accretion velocity enhances. 

The importance of magnetic diffusivity has been studied 
in several accreting systems, such as the protostellar discs 
(Stone et al. 2000; Fleming & Stone 2003), discs in dwarf 
nova systems (Gammie & Menou 1998), the discs around 
black holes (Kudoh & Kaburaki 1996), and Galactic centre 
(Meha & Kowalenkov 2001; Kaburaki et al. 2010). Moreover, 
two and three-dimensional MHD simulations have shown 
that resistive dissipation is one of the crucial processes that 
determines the saturation amplitude of the magnetorota- 
tional instability (MRI). As, the linear growth rate of MRI 
can be reduced significantly due to the suppression by ohmic 
dissipation (Sano et al. 1998; Fleming et al. 2000; Masada 
& Sano 2008). Moreover, from comparison of ideal and re- 
sistive MHD simulations, it seems the magnetic diffusivity 
may play an important role in astrophysical outflows (Fendt 
& Cemeljic 2002; Cemeljic et al. 2008). 

As mentioned, semi-analytical studies of ADAFs with 
outflow are typically related to systems without magnetic 
diffusivity. While, non-ideal MHD simulations imply that 
the resistivity plays an important role on outflows (e. g. 
Fendt & Cemeljic 2002; Cemeljic et al. 2008). Akizuki & 
Fukue (2006) proposed a self-similar solution for ADAFs 
with a highly ionized gas. Thus, they assumed that the 
plasma resistivity is zero, and only viscosity is due to tur- 
bulence and dissipation in the disc. Also, they ignored the 
effects of outflow/wind in their model. In this paper, we 
want to explore the effects of the resistivity on magnetized 
ADAFs in the presence of outfiows. Thus, we adopt the pre- 
sented solutions by Knigge (1999), Akizuki & Fukue (2006), 
and Shadmehri (2008). The paper is organized as follow. In 
section 2, the basic equations of constructing a model for 
ADAF in the presence of toroidal magnetic field, resistiv- 



ity and outfiows will be defined. In section 3, self-similar 
method for solving equations, which govern the behavior of 
the accreting gas will be used. Results of the present model 
are brought in sections 4 and 5. The summary of the model 
will appear in section 6. 



2 BASIC EQUATIONS 

We consider a non-ideal magnetohydrodynamics of steady, 
axisymmetric, viscous, accreting and rotating fiow in pres- 
ence of a purely toroidal magnetic field. We use a cylin- 
drical coordinate (r, 0, z) centred on the accreting object. 
We ignore the general-relativistic effects and use Newtonian 
gravity. Under these assumptions, the continuity equation is 



0, 



(1) 



where Vr is the radial velocity, E = 2pH is the surface den- 
sity of the disc, p and H being the mid-plane density and 
half-thickness of the disc, respectively, and Af^, is the mass 
loss rate by outflow/ wind. The half-thickness of the disc is 
given hy H = VI + H/fix, where Qk is the Keplerian 
angular velocity, Cs is the sound speed, H is deflned below 
in equation (9). It will be reduced to its traditional form of 
H — Cs/^K in absence of the toroidal component of mag- 
netic held {Btp-, equation 8). The sound speed is defined as 
Cs = (Pgas/ p)^^^ , where pgas is the gas pressure. The cumu- 
lative mass-loss rate from the disc can be written as 



M„(r-) 



4Tvr' mu,(r') dr' , 



(2) 



where ri„ denotes the radius at the inner edge of the disc 
and rhw{r) is the mass-loss rate per unit area from each disc 
face. Since the mass accretion rate is M = —2TTr'Evr, from 
the equations (1) and (2), we can write 



dM 
dr 



dr 



(3) 



Above equation implies that the mass accretion rate varies 
by radius due to outfiow. Thus, we exploit a power-law de- 
pendence for mass accretion rate as follows (e.g. Blandford 
& Begelman 1999) 



M{r) = M{R)(^' 



(4) 



where R is the radius at the outer edge of the disc, M{R) 
is the mass accretion rate at 7?, and s is a free parameter, 
which for a disc without outflow/wind, s = and in the 
presence of the outfiow/wind, s > (e.g. Fukue 2004). The 
observed broad-band spectra of Sgr A* and soft X-ray tran- 
sients can also be fitted by ADAF models with moderate 
outflows, s ~ 0.3 - 0.4, if the direct heating of electrons in 
ADAFs is efficient (Quataert & Narayan 1999; Yuan et al. 
2003). Equations (l)-(4) imply that 



rh^{r) 



M{R) 



47ri?2 \RJ 
The radial equation of momentum is 



(5) 



Vr—r- = r (Q^ - 
dr 



where fl is the angular velocity of the fiow and ca is Alfven 
speed, which is defined as Ca = B^/{4np) = 2p,nag/p, Pmag 
being the magnetic pressure. 
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The angular momentum transfer equation with consid- 
eration of the outflow/wind can be written as (e.g. Knigge 
1999) 



r dr 



r vT.— I - 
dr 



{Irfn dM^ 
27rr dr 



(7) 



where the two terms on the right-hand side of above equa- 
tion describe the effects of viscous torques due to shear {i/, 
which is the effective kinematic viscosity) and the angular 
momentum sink provided by the outflow. Here, it will be as- 
sumed that matter outflowed at radius r on the disc carries 
away speciflc angular momentum {lr)^Q,. Thus, I = corre- 
sponds to a non-rotating disc wind and / = 1 to outflowing 
material that carries away the specific angular momentum 
it had at the point of outflow (Knigge 1999). 

The hydrostatic balance in the vertical direction is in- 
tegrated to 



GM 2 2 



l + i(^ 

2\cs 



= {l + U)c 



Here, we introduce the parameter 11 by 



n : 



Pgas 



(8) 



(9) 



which is the degree of magnetic pressure to the gas pressure. 
Since we will apply a steady self-similar method to solve 
system equation, this parameter will be constant through- 
out the disc. Really, this parameter is a function of position 
and time (Machida et al. 2006; Oda et al. 2007; Khesali & 
Faghei 2008, 2009). Studies of hot accretion flows represent 
the typical value of H lies in the range 0.01-1 (De Villiers et 
al. 2003; Beckwith et al. 2008), but here we also consider the 
magnetically dominated case (H > 1). Because, MHD simu- 
lation by Machida et al. (2006) shows as thermal instability 
grows in an accretion flow, the magnetic pressure exceeds 
the gas pressure due to the disc shrink in the vertical di- 
rection and conservation of the toroidal magnetic flux. This 
will result in large H and forms a magnetically dominated 
accretion flow (Oda et al. 2007). 

We assume both of the viscosity and the diffusivity are 
due to turbulence in the disc, so that it is reasonable to use 
these parameters in analogy to the a-prescription of Shakura 
& Sunyaev (1973) for the turbulent 



Pm.ri = aCsH, 



(10) 



where Pm is the magnetic Prandtl number of the turbulence 
assumed a constant of order of unity, rj is the magnetic dif- 
fusivity, and a is a free parameter less than unity. 

Here, we can write the energy equation considering en- 
ergy balance in the system. We will assume the energy re- 
leased due to viscous and resistive dissipations can be bal- 
anced by the advection cooling and energy loss of outflow 
(e.g. Shadmehri 2008). Thus, 



'EVr dCg 

7 — 1 dr 



dr 



+ 



H 



27rP™ 



r,rr 2 dp 

dr 



r dr 



-2^m.uj{r)vjf{r), 



(11) 



where / is a constant less than unity and is called advec- 
tion parameter. The parameter / measures how much the 
flow is advection-dominated (Narayan & Yi 1994). The first 
two terms on the right-hand side of above equation repre- 
sent the energy generated due to viscous and resistive dis- 
sipation, respectively. The resistive dissipation is derived by 
{4tv / c'^)rj3^ , where J = (c/47r)V x B is the current density. 
Moreover, the last term on the right-hand side of energy 
equation is the energy loss due to wind or outflow. Depend- 
ing on the energy loss mechanism, dimensionless parameter ^ 
may change. We consider it as a free parameter of our model 
so that larger ^ corresponds to more energy extraction from 
the disc due to the outflows (Knigge 1999; Shadmehri 2008). 

The creation/escape rate of the magnetic field can be 
described by dynamo and diffusion. We define the advection 
rate of the toroidal magnetic field as (Oda et al. 2007) 



(12) 



which is used instead of the induction equation. Since we 
study a steady-state accreting system, the above quantity 
will be constant in absence of the dynamo and the diffusion 
effects. This quantity can vary with radius due to the pres- 
ence of the dynamo/diffusion effect (Machida et al. 2006; 
Oda et al. 2007). In the present model, we expect the mag- 
netic fiux advection rate varies with radius due to the pres- 
ence of resistivity. We will consider this property in next 
section. 



3 SELF-SIMILAR SOLUTIONS 

We seek self-similar solutions in the following form (e.g. Ak- 
izuki & Fukue 2006; Shadmehri 2008): 



E(r) = coEo 



Vr{r) 



GM 



C3 



fGM\ 



r \ ' 



2 / s 



2nc£ 



fGM\ 
\ R ) \R 



(13) 
(14) 
(15) 
(16) 
(17) 



where So and R are exploited to write the equations in 
the non-dimensional forms. Substituting the above solutions 
in the continuity, radial momentum, angular momentum, 
hydrostatic, and energy equations [(l),(6)-(8), and (11)], we 
can obtain the following relations: 



cqCi — m 



1 2 



C3 - + 1 = 



(18) 



(19) 
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ci - 3ac3(s+-)Vl + n 



+ / srh = 



H 



\/(i + n)c3 



(20) 



(21) 



Q!n/(7 - 1) ( S - ^ ) C0C3 + PmCoCs X 



9 ^,2 2ci 

2 ^ ' %/TTrT 



2smgP^(7- 1 

^/^TrI 



(22) 



where m is non-dimensional mass accretion rate and is de- 
fined 



MiR) 



2-kT.oVGMR 



(23) 



Using the equations (18)- (22), we obtain a quadratic 
equation for the coefficient of C3: 



We cannot extend the solutions for larger values of He, be- 
cause the right-hand side of equation (27) becomes neg- 
ative and a negative C2 is clearly unphysical. Moreover, 
He = 19>Pm implies that the flow can be magnetically dom- 
inated for Pm > 1/18 ~ 0.056. It means the flow can be 
had a strong magnetic field in the presence of resistivity. 
This property is in accord with resistive MHD simulations 
of Machida et al. (2006). 

Now, we can investigate the radial dependence of the 
magnetic fiux advection rate (equation 12). The self-similar 
solution for this quantity implies that 



'i'(r) 



^out (-^) 



(s-3/2)/2 



(29) 



where ^out is the magnetic flux advection rate at outer 
edge of the disc, R. Since s < 3/2, the magnetic fiux in- 
creases with approaching to central object and the stronger 
wind/outfiow reduces this increasing. The radial dependence 
of $ is qualitatively consistent with results of Machida et al. 
(2006) and Oda et al. (2007). 



1 + 2s 



2Ps - 1 



(7-1)/ 



1 + 2s 



2/2 s - 1 



+ 



9 
2 



1 + 



2s^ { l-f 2s 



3/ 

and the rest of the coefficients are 



Co 



Cl 



3 aVTTrT 



2rs - 1 



l + 2s 



C3 



0, (24) 



(25) 



2 

C2 




(26) 



(27) 



Without mass outflows, resistivity, and toroidal magnetic 
field, i.e. s = / = ^ = 0, Pm = 00, and 11 = 0, above rela- 
tions reduce to the standard ADAF solutions (Narayan & Yi 
1994). Moreover, in the absence of wind and resistivity but 
with toroidal magnetic field, above relations reduce to re- 
sult of Akizuki & Fukue (2006). However the present model 
includes toroidal magnetic field, outfiows, and resistivity. 

The studies of resistive and magnetized ADAFs (e. g. 
Faghei 2011) imply that the solution for a set of the input 
parameters reaches to a non-rotating limit at a specific of 11 
which we call it by lie. Assuming C2 = and s = / = ^ = 
(no wind case) in equations of (18)-(22), lie can be written 
as 



18P„ 

/ 



5/3-7 



(28) 



For typical values of adiabatic index and advection parame- 
ter in ADAFs, 7 = 4/3 and / = 1, we can write He = 18Pm. 



4 RESULTS 

Now we can investigate the behavior of the solutions in the 
presence of the outfiows and resistivity. The effects of out- 
fiows and resistivity are studied via parameters s, I, ^, and 
P^^ ■ Here, the inverse of Prandtl number specifies the re- 
sistivity of the fluid. Because, 77 oc a/Pm and a parameter 
is 0.1 in all Figures. The behavior of physical variables as 
a function of P^^ are shown in Figures 1 and 2. The solu- 
tions in Figures 1 and 2 represent the radial inflow speed, 
Cl, and sound speed, C3, both increase with the magnitude 
of resistivity (P,~^). These properties are qualitatively con- 
sistent with Faghei (2011). The density profiles, cq, show 
that it decreases by adding resistivity. It can be due to tem- 
perature raise of the flow. The rotational velocity decreases, 
C2, with the magnitude of resistivity. Because, the viscous 
torque increases with the temperature (;/ oc oc T). These 
properties are also consistent with results of Faghei (2011). 

In Figure 1, we also studied the effect of parameter s 
on physical variables. The value of s measures the strength 
of outflow and a larger s denotes a stronger outflow. Fig- 
ure 1 shows that for non-zero s, surface density is lower 
than the standard ADAF solution and for stronger outflows 
this reduction of surface density is more evidence. We can 
see that ADAF with wind rotates more quickly than those 
without winds and lead to enhance accretion velocity. The 
solution shows that temperature decreases for stronger out- 
flows. On the other hand, outflows play as a cooling agent. 
These properties are in accord with results of Shadmehri 
(2008) and Abbassi et al. (2010). 

In Figure 2, the effect of energy loss due to outflows 
is studied by ^ parameter. As with the magnitude of ^ pa- 
rameter, the more energy will carry by outflows. Due to this 
energy loss, we expect the temperature should decrease by 
adding the ^ parameter. Temperature profiles confirmed this 
property. Since the turbulence viscosity is proportional to 
temperature (iv oc T), the efficiency of angular momentum 
transport decreases with the ^ parameter. Decreasing vis- 
cous torque increases the rotational velocity and decreases 
the radial infall velocity. These properties are consistent 
with previous works (e. g. Shadmehri 2008). 
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Figure 1. Physical quantities of the flow as a function of Pm 
for 7 = 4/3, a = 0.1, H = 0.5, and / = « = g = 1. The sohd, 
dotted, dashed, and dot-dashed hues represent s = 0, 0.01, 0.02, 
and 0.03. 



The physical variables as a function of parameter 11 
and several values of mass and energy losses are shown in 
Figures 3 and 4. By adding 11 which indicates the role of 
magnetic filed on the dynamics of accretion discs, we see 
the sound speed becomes larger, while surface density de- 
creases. Moreover, Figure 3 represents that the radial and 
rotational velocities increase with the magnitude of H. In- 
crease in the radial velocity is due to the magnetic tension 
term dominates the magnetic pressure term in the radial 
momentum equation, which assists the radial velocity of ac- 
cretion flows. Moreover, increase in the rotational velocity 
is because of that the disc should rotate faster than the case 
without the magnetic field which results the magnetic ten- 
sion. These properties are qualitatively consistent with the 
previous works on magnetized ADAFs (e.g. Akizuki & Fukue 
2006; KhesaU & Faghei 2009; Abbassi et al. 2010). Figures 3 
and 4 imply that the mass and energy losses due to wind give 
the same results by Figures 1 and 2. On the other hand, the 
effects of mass and energy losses on the physical variables 
do not change in low and high values of magnetic field. 

Figures 5 and 6 represent that the disc thickness in- 
creases with the magnitude of the resistivity or the magnetic 
field. Because the temperature increases by adding the resis- 
tivity or the magnetic field, and from equation (21) we can 
see the disc thickness increases with temperature. In Fig- 
ures 5 and 6, the disc thickness is also studied for several 
values of parameters a and The disc thickness profiles im- 
ply that it becomes thinner for the stronger mass or energy 
losses due to outflows. Because these parameters reduces the 
temperature of the flow. 




Figure 2. Same as Figure 1, but s = 0.1, and the solid, dotted, 
dashed, and dot-dashed lines represent § = 0,0.1,0.2, and 0.3. 



5 THE BERNOULLI PARAMETER 

Here, we exploit Bernoulli parameter (Be) to consider the ef- 
fects of the resistivity to generate/enhance outflows in mag- 
netized ADAFs. Because, this parameter measures the like- 
lihood that outflow or wind may originate spontaneously 
(Narayan & Yi 1994). An adiabatic flow has a constant Be 
along streamlines. If Be is positive for any of accreting gas, 
then this gas can potentially reach infinity with a net posi- 
tive kinetic energy. The Bernoulli parameter, defined as the 
sum of the kinetic energy, the enthalpy and the potential 
energy of the accretion flow, 

hv^ +^n^) + ^cl~r^nl. (30) 



Be ■- 



7- 



1 



Narayan & Yi (1994) showed that Bernoulli parameter 
is positive in height-integrated advection dominated flows, 
and suggested this may explain the frequency occurrence 
of outflows and wind in many accretion systems. Using the 
self-similar transformations of (13)-(17), the Bernoulli pa- 
rameter can be written as 



h = 



Be 



7 



2(ci+C2) + ^_^ 



C3 



(31) 



where h is the normalized Bernoulli parameter. In Figure 
7, the behavior of this parameter as a function of P,~^ is 
studied for different values of magnetic field. The Be profiles 
represent that it is positive and increases with the magnitude 
of resistivity or magnetic field. It can be due to increase of 
the flow temperature by adding the resistivity or magnetic 
field. Moreover, the Be profiles show the magnetic field is 
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Figure 3. Physical quantities of tiie flow as a function of 11 for 
7 = 4/3, a = 0.1, P„ = 1, / = 1, / = 1 and g = 1. The solid, 
dotted, dashed, and dot-dashed linos represent s = 0, 0.01, 0.03, 
and 0.05. 



more important in high resistivity. Thus, the outflows can 
be enhanced in resistive and magnetized ADAFs. 



6 SUMMARY AND DISCUSSION 

Mass loss appears to be a common phenomenon among ac- 
cretion systems. The observational evidences of ADAFs im- 
ply that outflow is important in such system. Moreover, the 
non-ideal MHD simulation results represent that resistivity 
can play an importance role in outflow of accretion discs. 

In this paper, the structure of a magnetized ADAF in 
the presence of resistivity and outflow is investigated. We 
assumed that the magnetic fleld has a purely toroidal com- 
ponent. The outflow emanating affects on the equations of 
continuity, angular momentum and energy, and can there- 
fore act as a sink for mass, angular momentum and energy. 
We adopted the presented solutions by Knigge (1999), Ak- 
izuki & Fukue (2006), and Shadmehri (2008). Thus, we as- 
sumed that angular momentum transport is due to viscous 
turbulence and the a-prescription is used for the kinematic 
coefficient of viscosity. We also assumed that the flow does 
not have a good cooling efficiency and so a fraction of energy 
accretes along with matter onto the central object. To solve 
the equations that govern the structure behavior of magne- 
tized ADAF with outflow, we have used a steady self-similar 
solution. 

The present model represented the outflows in accre- 
tion flows can improve by the resistivity and magnetic field. 
These properties are in accord with resistive MHD simula- 
tions of Fendt & Cemeljic (2002) and Cemeljic et al. (2008). 



Figure 4. Same as Figure 3, but s = 0.1, and the solid, dotted, 
dashed, and dot-dashed hues represent g = 0, 1, 2, and 3. 
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Figure 5. The ratio of height thickness to radius as a function 
of for different s (left panel), and ^ (right panel). The input 
parameters in left and right panels are same as Figures 1 and 2, 
respectively. 



For example, Fendt & Cemeljic (2002) showed that resistiv- 
ity can affect the outfiows structure in accretion discs. More- 
over, they found that the outfiow velocity increases with the 
magnitude of the resistivity and toroidal magnetic field. 

In this paper, we assumed a purely toroidal magnetic 
field that provides restrictions in the present model. For ex- 
ample, a purely toroidal magnetic field is not enough to have 
magnetically driven outfiows. Thus, the present model is un- 
suitable to use in disc with magnetically driven outfiows. 
Moreover, we considered the presented model in a height- 
integrated approach and applied it to calculate Bernoulli 
parameter. Narayan & Yi (1995) showed that Bernoulli pa- 
rameter in ADAF varies by latitude. As Bernoulli parameter 
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Figure 6. The ratio of lieight thickness to radius as a function 
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Figure 7. The dimcnsionless of Bernoulli constant as a function 
of for different values of 11. The input parameters arc same 
as Figure 1, but s = ^ = ( = 0, and the solid, dotted, dashed, and 
dot-dashed lines represent 11 = 0.1, 0.3, 0.5, and 0.7, respectively. 
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in equatorial is negative and becomes positive for high lati- 
tude. One can investigate latitudinal behavior of the present 
model. 
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